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The "equations of motion" of a linear electroacoustic transducer arc written in the alternative 
forms A\Y = 0, /, ,Y = 0, where the matrices L 1 are linear differential expressions, X = {E t //, u) rep- 
resents electroacoustic fields, and the superscripts distinguish selected normalizations of the equa- 
tions of motion. To each operator L- corresponds a mathematically defined adjoint operator {L' )" 
and an associated adjoint transducer, whose material tensor parameters arc given by certain trans- 
positions and interchanges of the parameters of the given transducer. Dissipative characteristics 
(lossiness, losslessness, or "gaininess") of the material of the given transducer arc preserved point- 
wise in the adjoint transducers. A «£eneralized reciprocity lemma leads to relations of reciprocal type 
between external properties of the given and the adjoint transducers. In the self adjoint cases, the 
conventional electroacoustic reciprocity and antireciprocity relations are obtained and the derivation 
of those relations is critically confirmed. The generalized or adjoint reciprocity relations have been 
applied in the plane-wave scattering-matrix formalism developed for electroacoustic transducers bv 
D. M. Kerns. 

Corollaries of the adjoint reciprocity relations, analogous to conventional reciprocity theorems, 
but involving properties of adjoint pairs of transducers, are derived. Examples are discussed of trans- 
ducers for which it is feasible to form the adjoint transducers. 
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Glossary of Symbols 

Denotes the adjoint. 

Denotes the transpose. 

Surfaee of transducer; surface of power supply or detector; and surface of trans- 

ducer-plus-(power supply or detector). 

Amplitude of incident and emergent em mode which feeds the transducer. 

Incoming and outgoing plane-wave amplitudes in ambient fluid. 

Direct and converse piezoelectric stress tensor (triadic). 

Amplitude of spherical wave in ambient fluid. 

Harmonic and static magnetic induction. 

Direct and converse piezomagnetic stress tensor (triadic). 

Hooke's tensor (tetradic). 

Electric displacement. 

Harmonic and static electric field. 

Electric and magnetic field components of em mode which feeds the transducer. 

A MS Subject Classification: 7000. 7800. 
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e: Permittivity dyadic. 

r) = yla)p : Plane-wave admittance. 

170 : Admittance of em mode feeding transducer. 

F: Bilinear concomitant associated with the operator L. 

f,f s : Induced and applied volume forces. 

y: z-component of the plane-wave propagation vector. 

H : Magnetic field. 

h, h'\ Spacial impedances describing transducer. 

/: Identity dyadic. 

7, J s : Induced and applied current densities. 

K, L: Transverse part of kand /. 

A\ /: Plane-wave propagation vectors in ambient fluid. 

L, L~: Transducer operator matrices. 

M s : Applied magnetization. 

M(f): Microphone response. 

p: Permeability dyadic. 

p<>: Free-space permeability. 

p, p () : Harmonic and static pressure. 

P s : Applied polarization. 

p<>: Static mass density of the ambient fluid. 

p m , p° M p' m \ Harmonic, static, and total mass density of transducer material. 

p r , p° : Harmonic and static charge density. 

So: Transducer feed surface. 

Si: Reference xy plane in ambient fluid. 

Soo, S01, S10, Sn: Plane-wave scattering matrices describing transducer. 

S(r): Speaker response. 

5: Strain dyadic. 

cr : Conductivity dyadic. 

T 9 To, T': Harmonic, static, and total stress dyadic. 

t, v\ Magneto-electric dyadics. 

u: Fluid or material velocity. 

Vo, /o: Input voltage and current to transducer. 

V{K), I{K): Voltage and current spectral densities. 

Zft, Zo : Spacial impedances (along with h, h' ) describing transducer. 

1. Introduction 

Analytical techniques developed by Kerns [1, 2, 3] 1 for the measurement of microwave antennas 
have been "translated" recently by that author into corresponding techniques for -lie measure- 
ment of electroacoustic transducers in fluids. The basic theory is formulated with a scattering- 
matrix description and emphasizes the use of plane-wave spectra for the representation of 
acoustic fields. For a summary and for details of the theory and its variety of applications, reference 
should be made to Kerns [4]. 

Although the electromagnetic (em) theory is quite comprehensive, the acoustics applications 
are strengthened by generalized or adjoint reciprocity relations [5], the derivation of which forms 
the principal subject of the present paper. The generalized reciprocity relations are expressed in 
the notation of the Kerns plane-wave scattering-matrix (PWSM) description, and also in the 
notation of the spacial impedance-matrix (SIM) description introduced by Foldy and PrimakofT 
[6, 7] in their classic papers on electroacoustic transducers and reciprocity. The PWSM expres- 



1 Figures in brackets indicate the literature references at the end of this paper. 
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sions of adjoint reciprocity have become especially important because of their use in the theo- 
retical foundation of extrapolation techniques which predict gain on an arbitrary axis from measure- 
ments at reduced distances along that axis [4]. (Extrapolation techniques were introduced for 
antennas by Wacker and Bowman [8]. ) 

The following approach is taken to establish the adjoint reciprocity relationships. The equa- 
tions of motion and constitutive relations describing a large class of linear electroacoustic trans- 
ducers are reduced to the form of a linear operator matrix. The operator is "transposed" to form a 
unique adjoint operator, which is compared and combined with the original operator to derive a 
generalized electroacoustic reciprocity lemma relating the fields of the given transducer to those of 
an adjoint transducer. The material parameters of the adjoint transducer are determined by certain 
transpositions and interchanges of the material parameters of the given transducer. The reciprocity 
lemma leads to the generalized reciprocal relations between the characteristic (PWSM and SIM) 
matrices of the given and adjoint transducers. In addition, corollaries of the adjoint reciprocity 
relations analogous to conventional reciprocity theorems are derived for adjoint pairs of transducers. 

Power relations and the associated "value classes" are investigated also for electroacoustic 
transducers and their adjoints. In particular, it is shown that the dissipative characteristics (lossi- 
ness, losslessness, or "gaininess") of the given transducer arc preserved pointwise in the adjoint 
transducers. 

Although a unique adjoint corresponds to a given operator, many different operators can rep- 
resent the equations of the same transducer. For example, multiplication of one of the equations 
of motion by a constant will change the operator matrix without changing the content of the equa- 
tion of motion. Analysis shows, however, that for a given transducer only two operators, L + and 
L (from here on the "-h" will be omitted from L* ), have an associated adjoint operator, L n and A", 
which represents the equations of another transducer (the adjoint transducer). A transducer self- 
adjoint (L" = L or L "=L ) with respect to L or L corresponds to the conventional reciprocal or 
antireciprocal transducer respectively. One of the advantages of the linear operator approach is 
that it facilitates the identification of self-adjoint transducers. 

A number of experimental measurement techniques applicable to reciprocal (or antirecipro- 
cal) transducers only, could be extended to nonreciprocal transducers if the adjoint transducers 
were producible in the laboratory. In general, the difficulties in synthesizing the necessary electro- 
acoustic material parameters may discourage or prevent the creation of the adjoint transducers. 
However, two groups of transducers are discussed for which the creation of an adjoint transducer 
involves, in principle, a simple modification of the given transducer. 



2. The Transducer System 

The electroacoustic transducer system under consideration is pictured in figure 1. The trans- 
ducer is bounded by the closed, finite surface Aq. The em power supply or detector, which is 
bounded by the closed, finite surface A f) , feeds the transducer through the area So common to both 
A{) and A,,. The transducers, which may contain static bias fields and their sources, are termed 
"passive" if they cannot radiate more power than they absorb. The present analysis, however, is 
not limited to passive transducers, but may include transducers which are active because they con- 
tain "gainy" material, or radiating sources, or both. (A transducer will be called source-free if its 
internal equations of motion are homogeneous, i.e., if they are satisfied by all fields equal to zero.) 

The transdueer-plus-( power supply or detector) is immersed in a homogeneous, isotropic, sta- 
tionary fluid which extends to infinity and supports the time harmonic (e il °', co > 0) acoustic 
pressure-velocity field of small (first order differential) amplitude. The static pressure, static den- 
sity, and propagation constant (the latter two of which may be complex function of real frequency co 
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to account for viscous and "expansive friction" losses) are constants in space and time denoted by 
Po, po and k respectively. The fluid between the infinite xy plane S\ (z= ) and the external surface 




Po^o.k 



FIGURE 1. Schematic of the transducer system. 

A=A„ + A P denotes the external surface of the electroacoustic transducer-plus! power supply or detector). 



a(K) 



-► b(K) 



A of the transducer-plus-(power supply or detector) contains no acoustic sources, 2 although arbi- 
trary sources (at frequency a> ) may exist to the right of Si. 

The normal velocity on the surface A v of the power supply-detector, including the feed area 
So, is assumed negligible. Electromagnetic shielding necessary for deriving the adjoint reciprocity 
relations is discussed in section 5. 

The area So may designate the perpendicular cross section of an open or closed em waveguide 
(uniform and isotropic) in which a single em mode is propagating with incident amplitude a<> 
and emergent amplitude 6o, or simply the cross section of two wire leads when the frequency a> 
is low enough to allow analysis by quasistatic voltages and currents. The cross-sectional area So 
for a mode on an open waveguide (such as a Lecher line or microstrip) must extend sufficiently 
far from the guide to insure that the fields in the mode outside So contribute negligibly to the 
normalization integral ((7) below). The surfaces A p , So, and A {) represent imaginary boundaries 
in that they art 1 chosen to facilitate the theoretical analysis and need not coincide with the 
physical boundaries of the system. Of course, the surface A— Aq-\- A p must be chosen within 
(or on the interior boundary of) the ambient fluid. 

The continuity and momentum equations, which govern the harmonic acoustic fields in 
the source-free fluid between the surface A and the infinite plane Si, may be expressed as 



k 2 p = —i(t)poV - u 
Vp = icop o/7 , 



(la) 
(lb) 



2 BeU 
i. b\ a 



A and S, there 
.metric Hooke's 



an he finite regions of ei 
ensor. However, the rele 



r acoustically 
ce plane S, n 



terial or m< 
lin entirely 



rial in which volume forces are zero 
the homogeneous ambient fluid. 



id stress is linearly related to 



20 



where the complex amplitudes, p and /7, denote the position-dependent excess pressure and 
fluid velocity. The effect of gravity on harmonic variation in mass density is neglected. 

The excess pressure and normal velocity on the surface Si may be expanded in a double 
Fourier integral (transform) over plane wave voltage and current spectral densities [4]: 

P(fi)-±f V(Z)e*- M diC (2a) 



*<*)—£/ 



I(K)e lKR dK. (2b) 



Conversely, the plane-wave voltage and current densities, which are defined in terms of the 
incoming a(K) and outgoing b(K) acoustic plane-wave amplitudes, 

V(K) =a(K)+b(K) (3a) 

I(K) = [a(K)-b(K)]r)(K), (3b) 

may be expanded in a double Fourier integral (inverse transform) over the excess pressure and 
normal velocity: 

V(K)=-^f s P(R)e~** dR (4a) 

l(K)=-^f uAR)e~' R ■* dR. (4b) 

J S { 

The vector K = K. r i + K i J can be interpreted as the transverse part of the plane-wave propagation 
vector k. The vector/? is confined to the xy plane Si, and r)(K) equals y/cupo, where y is defined as 
(k 2 — K 2 ) 1/2 with the sign of the radical chosen to keep the real and imaginary parts of y positive. 

The transverse em fields (£,. H, ) on the cross section So of a waveguide may be written as 
the sum of the transverse electric e and magnetic ho components of the mode traveling to the 
right with incident amplitude a<) and to the left with emergent amplitude bo: 

E t = Voio, H t = Ioho, (S) 

where 

V { ) = a {) +bu (6a) 

/«>= (ao — Mrjo. (6b) 

If the dimensions of Co and h {) arc chosen as (meter) -1 and (eo, rjoAo) are consistent with 
Maxwell's equations, Vo behaves dimensionally as volts, / as amperes, and rjo as an admittance 
(ohm -1 ) normally chosen as a positive real number. Moreover, normalization may be expressed 
as a nondimensional number equal to unity, i.e., 

I (e () X/i () ) ■ £da = l, (7) 

is'o 

where h is the inward normal to the transducer. (The rationalized MKS system of units is used 
throughout.) If the area So simply cuts two wire leads at quasi-static frequencies, Vo and /o refer 
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to conventional circuit voltages and currents which do not necessarily serve as genuine modal 
coefficients. In that case, (6) become a definition of ao and bo. 

The PWSM description relates the outgoing amplitudes, bo and b{K) , to the incoming am- 
plitudes, a<» and a (K) , by the scattering equations (see ref. [4]), 

bo = Sooa« + I S lH (L)a(L)dL (8a) 

b(K)=Sio(K)ao + f Su (Kl)a(L) dl (8b) 

The scattering matrix is resolved into the four submatrices Soi, Sio, Su and Soo which represent 
the receiving, radiating, acoustic scattering, and em reflection properties respectively for the 
electroacoustic transducer. Like K, the vector L can be interpreted as the transverse part of the 
propagation vector/. 

Alternatively, the SIM relations 3 characterize the external behavior of the transducer by 
relating the "voltage" Vo and excess pressure /; on the surface of the transducer A {) (except possibly 
on the feed area So) to the "current"' /<» and the inward (to transducer) normal velocity u n on the 
surface Ao: 

V{) = Zj/o-T- h f (f )tt„(fo) da {) (9a) 

Ja 

p{r) = h(r)I {) + Zo(r,r )u n (fo) dao. (9b) 

Ja o 

Although three of the four spatial impedances (Z&, h' , h, Z ) depend upon the position vectors 
(r,f {) ) to the surface A () and all four depend upon the harmonic time frequency, they are assumed 
independent of the medium or sources surrounding the transducer and, as linearity implies, 
independent of Fo, /o,p and u„. 

It can be shown [9], without reference to the internal behavior of the transducers, that the 
PWSM equations (8) imply the SIM relations (9) but not vice versa. That is, the PWSM description 
applies to a larger class of transducers than do the SIM relations. However, the alternate ex- 
pressions of reciprocity (and antireciprocity), 

r)oSoi(K) = ± i?(/OSio(- K) (10a) 

r 1 (K)Sn(K,L)=r } (L)Su(-L,K), (10b) 



h'(r)=±h(r) (Ha) 

Z (f, fo)=Z (f , r) [6, 7], (lib) 

which are derived in section 5, can be proven [9] equivalent for transducers which obey both 
descriptions. 

In analogy to the 2-port cases discussed by McMillan [10], Kerns [4] defines reciprocity, anti- 
reciprocity, and nonreciprocity for electroacoustic transducers in terms of whether (10a) holds 
with the plus sign, the minus sign, or not at all. Foldy and'Primakoff [6,7] do not use the term 
"antireciprocity" explicitly, but distinguish between the plus and minus sign in (11a) by con- 



3 Foldy and Primakoff refer to equations (9) as the "linearity relations.*' Here they are referred to as the "spacial impedance-matrix relations'* to avoid confusior 
with the spectral scattering-matrix equations (8). which are also linear. 
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sidering reciprocity through electric-type coupling and magnetic-type coupling. 

Equations which characterize the internal behavior of the eleetroacoustic transducers under 
consideration are introduced in the following section. 

3. The Operator Matrix Description 

Primakoff and Foldy [7], in their original work on eleetroacoustic reciprocity, examined the 
equations which model the internal operation of a large class of linear transducers. They found 
that if certain "symmetry conditions" were satisfied by its material parameters and static fields, 
the transducers obeyed the SIM reciprocity relations (11). The present derivation begins in a 
similar way by writing the linear, harmonic (e - lojt ) equations that govern the "physics" of the mate- 
rial within the transducer A {) : 

V •f+i«p o w fi+/=0 (12) 
VxE-iwB = (13a) 

V xH+ia)D-J=0. (13b) 

The momentum equation (12) expresses Newton's second law for stationary (/7 = 0) media, and 
(13) express Maxwell's equations in differential form for stationary media in which regions of 
electric and magnetic polarization may exist. 

The static quantities, which will be distinguished by the superscript or subscript "0," are 
real functions of the space coordinates. All other quantities, except the real frequency to, represent 
complex, space-dependent, harmonic amplitudes, which superimpose upon the static variables as 
products with e~ fftrf . For example, if the total stress dyadic at a point in space and time is denoted by 
T' (r, t) then T' can be divided into a real static part T {) (T) and a complex harmonic part T(T)e~ la)t 
i.e., 

f{r,t) = f«(f) + f(f)e- ,wr . 

Of course, it is the real part of T' that corresponds to a measurable quantity. Only the complex 
harmonic amplitude T(r) shows in (12). Similarly the total mass density p' m (r, t) may be divided as 

Pm(r 9 t)= p&(f) + p w (f)e-' w/ , 

but only the real static part occurs in (12). In fact, the mass continuity equation is not included 
in addition to the momentum equation (12) because it merely introduces the extra variable p m (r). 

The other variables in (12) are zZ, the harmonic velocity of the material, and/ the total har- 
monic volume force. The usual harmonic em field vectors are denoted by £, //, #, Z), and J. All 
variables refer to macroscopic quantities which are evaluated, at least in concept, by averaging 
microscopic variables over small but finite volumes and time intervals. 

The number of unknowns in (12) and (13) reduces to the number of equations if the following 
linear constitutive relations characterize the transducer material: 

D= e'l + rff + o + P s (14) 

B= v-E + pH + $'3+ inM H (15) 

J=a-{E + a x/? ) + p°efi + Js (16) 

T=c--E'E-JI'H (17) 
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f=p° e E + JxBo+f s , (18) 

where the strain dyadic J is related to the velocity by 

s = -}-fou-l (VXu) X/V (19) 

io) \ Z J 

(I is the identity dyadic). 

The various "dot products" involved in the vector and dyadic transformations of (14-19) can 
be understood in terms of Cartesian tensor notation, in which €•£, 5' • E, a : 1 and W:f become 
eijE j, a'jjkE/-. ctjjA-SjA-, an d CtjkiSki respectively (summation over repeated dummy indices is im- 
plied). The "/z-adic" notation [11] is used wherever convenient because it preserves the familiar 
vector-dyadic notation for the electroacoustic fields [7, 12]. 

The dyadics o\ €, /u, and (f, V) stand for the conductivity, permittivity, permeability, and 
magneto-electric properties [13] of the transducer material. 

The direct piezoelectric and piezomagnetic stress triadics [14], a and /3, may be chosen 
symmetric: 

otijk = ot ik j (20a) 



ptjk = Ptkj, (20b) 



since the strain dyadic is symmetric. 



— 1 f diij , du k 1 _„. 

2ioj { dx k dxj J 

Similarly, because the stress is a symmetric dyadic (Tjj=Tji) (provided there are no distribu- 
tions of body torque, stress couples, or angular momentum density [15] ), the converse piezoelectric 
and piezomagnetic stress triadics [14], a' and /3\ must possess the symmetry, 



l jik 



(22 a) 



P'tjk = Pji*' (22b) 

The symmetry of the stress dyadic also requires the Hooke's tetradic c to obey the symmetry 

ctjki= Cji k i, (23a) 

and the symmetry of the strain dyadic allows the defined symmetry 

Cij k i = Ciji k . (23b) 

The static charge distribution pj!, the static magnetic induction Bo, all other static quantities, and 
material parameters which may depend upon frequency a>, remain independent of the values of 
the harmonic variables. 

Although only transducers which obey the homogeneous equations are considered finally 
(source-free transducers), the harmonic sources of polarization P s , magnetization M s (/x is the 
permeability of free-space), current J s , and volume force f s have been included in (14), (15), (16), 
and (18) to allow insight into applications involving these sources. 

A number of other assumptions underlie the development of (14-19). In the current equation 
(16) the "Hall effect" of the magnetic induction on drift velocity is ignored, as well as harmonic 
variations in conductivity which might be found, for example, in a biased carbon microphone. 
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The volume forces on electric and magnetic dipoles, static currents, and harmonic variations 
in mass density (acted upon by gravity) have been omitted in the force equation ( 18). The effects 
of gravity are usually minute, and if the material that contains the dipoles or static currents is rigid 
or has insignificant effect on the essential behavior of the transducer, the volume forces on the di- 
poles and static currents are justifiably neglected also. The magnet in a dynamic speaker may 
experience appreciable volume forces, but because the magnet is relatively rigid, (12-19) describe 
such speakers. However, the operation of a hypothetical transducer containing dipole or static 
current elements which move through strongly inhomogeneous em fields would not, in general, 
be described by (18). 

Application of Maxwell's stress tensor reveals that forces on surface charges and currents can 
be evaluated by substituting the appropriate surface delta functions into (18) if nonconductors 
contain no static surface charge, and if the electric field just outside the surface of conductors, on 
which static surface charge exists, is used in the evaluation of p®E. The force term p e Eo is not 
needed in (18) because it vanishes in nonconductors where the harmonic charge distribution p e 
must be zero, and in conductors where the static electric field Eo must be zero. The standard con- 
denser microphone is well described by (12-19). 

The constitutive 1 equations (14-19) may be used to reduce the three physical laws (12), (13a), 
and (13b) to functions of three vectors, the electric field E, the magnetic field //, and the velocity 
field u. The final three equations are conveniently displayed in operator matrix form as 



{icoe - —a 



(V -X-icov-) 



(Vx-hicoT- ) (-a: V+ (aXB {) ) • -p5!) 



(-ico/Z-) (0: V; 



(V -a' -fBoXo 1 -p?) (V •/*' •) -(BoXaXBo) ■ + p e °BoX-i<opl 



10) 



II 



js-icoPs 
i0)jJLuM s 



(24) 



Since (VX/Z) X/is antisymmetric, its dot products with the symmetric 3, /3, and c vanish. All 
operator symbols of the 3X3 matrix containing three vectors are evaluated in succession from right 
to left. In full scalar matrix notation, the left side of (24) becomes a 9 X 9 matrix operating on a 9- 
element column matrix, and the right side another 9-element column matrix. 

The particular arrangement of elements in the matrix of (24) was chosen because it exhibits 
Maxwell's equations as a 2X2 sub matrix (the upper left corner), and because the associated 
adjoint operator matrix also represents the equations of an electroacoustic transducer (the ad- 
joint transducer). 



4. The Adjoint Operator Matrix and an Electroacoustic Reciprocity Lemma 

The development of the desired generalized reciprocity relations requires the intermediate 
derivation of a bilinear divergence expression which combines, essentially, the Lorentz lemma [16) 
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of electromagnetic theory with the acoustical reciprocity theorem ((45) below). The theory of 
adjoint operators assures the existence of such a bilinear expression (called the "bilinear iden- 
tity"), and provides a straightforward method for obtaining the identity from the operator matrix 
of (24) [17]. Specifically, if /. denotes the linear operator matrix of (24), the bilinear identity may 
be written in differential form as 

4>uL4>2-4>2tL a <l>i = V ' F> (25) 

where L" denotes the adjoint operator associated with L\ (f>\ and (/> 2 are any two continuously 
diffcrentiable column matrices containing three vectors (the subscript t transposes the column 
matrix to the corresponding row matrix); and V • F is the required divergence expression. The 
explicit evaluation of the bilinear concomitant F requires the determination of the adjoint op- 
erator L'\ which is actually defined by (25) [17 1. 

For a given linear operator L, the adjoint operator/." exists uniquely [9], and may be found by 
transposing the matrix elements of L, and replacing the differential operators by their adjoints. 4 
Fortunately, the adjoints of the four different linear differential operators of L, 

VX,S :V,V -5' -,andV -c : V, (26a,b,c,d) 

are readily determined as 

VX,-V-5r ,-"; :V,andV-* :V, (27a,b,e,d) 

respectively, with the help of the associated identities, 

A x • (Vxi 2 )-I 2 •(Vxi 1 ) = V- (A 2 xAi) (28a) 

A x • (V • {a -A,)) +A t ■ (5, : VA X ) = V - (A, -5-A 2 ) (28b) 

Ay • (V • (f:VA 2 )) -A 2 - (V- (S,±V/i)) = V- {A x -fiVAi-Az-ftiVAx). (28c) 

The subscript t denotes the transposed matrix, and ^1,^2 afe arbitrary, continuously differentiable 
vectors. In tensor notation the transpose of a ijk and c#jw, for example, used in (27), (28) and below 
is otjkj and Ckiij. Note that the curl operator is self-adjoint, that the direct and converse piezo- 
electric and piezomagnetic operators are adjoints of each other if at — — a', ft = — /3' r and that 
the purely acoustic operator is self-adjoint if Ct = c. These self-adjoint properties are understandable 
from Lanczos' theory of adjoint operators where he shows that the curl operator is self-adjoint 
and that the gradient and divergence operators are negative adjoints of each other [17]. 

The straightforward transposition of the nondifferential elements of L completes the con- 
struction of the adjoint operator L" : 

(iioec - av) (VX - ioyvr) -S/:V- cr*Xfi ) • - p\\ 

L«=f (Vx + icur,:) (-itojiti) (~|;:V) (29 ) 

(V-Sr-floXffr-p?) (— V-j5 r .) I \—^ (B {) x&,xB {) )--p {) B»x-ia) P l . 

J L 'W e J 

Comparison of the adjoint operator L a with L reveals that, at a given frequency oj, the adjoint 



4 Determination of adjoint differential operators as implied by the bilinear identity (25) is not synonymous with the interchanging of rows and columns of the 
Cartesian matrix form of the operator. Contrast, for example, Lanczos, ref. [17] with Auld, ref. [12]. 
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operator describes the equations of a second transducer (the adjoint transducer) in which the 
transposed parameters 



(T,, €r, flt.—Tt, 



replace the original parameters 



/3„";,-js;, 



B { 



,p; f andpj 



cr. e. /x, *>, r, 3', /3', a, /3. c, Z?<>, p|! and piJ, , 



(30a) 



(30h 



respectively. The adjoint transducer refers to a mathematically conceived transducer with hy- 
pothetical material parameters which may or may not be physically realizable (see sec. 7). When 
the transposed parameters (30a) equal the original parameters (30b), the operators L and L" be- 
come identical to form a self-adjoint operator for which the adjoint transducer and original trans- 
ducer are one and the same. In that case Bo is zero (/?<> = — Bo), which corresponds to the electric- 
type coupling of Primakoff and Foldy [7|. 

The expression (29) for the adjoint operator allows the explicit evaluation of V • F. If 02 
is restricted to solutions of (24), and $i to solutions of the adjoint equations corresponding to 
(24) with L" replacing/., i.e., 



L<j> 2 = Q,L"<l> l = Q n 



(33a, 



Q = 





E 




E" 




4,2 = 


H 
u 


. </>,= 


H" 

a" 




~js - ioP, 




~J" s -i(oP'' 


tw/xoM.s 


, Q" = 


iofitM", 


/« - So x J, 




f" + &« x J" 



(31a,b) 



(32a,b) 



(33a, b) 



then V • F becomes, after rearrangement and use of (17) and (19), 

V F = V ■ (EXH" -E"XH)-V ■ (u" f-u- f") 

= E" ■ (Js - io)P s ) + ioj/xM" ■ Ms + u" • (/, - Bo x Js) 
- E ■ {]'• - iwP") - iuixnH -M«-u- (f« + BoXj"). 



(34) 



When the harmonic sources are set equal to zero (in the transducer not the power supply), applica- 
tion of the divergence theorem to (34) in the volume bounded by the surface Ao of the transducer 
yields a generalized reciprocity ''lemma" between the independent electroacoustic fields of a 
given transducer and those of its mathematical adjoint: 



[ (EX H" 

'A 



I 



E" x H) ■ hda + (T" ■ a-T ■ a") ■ hda = 0. 



(35) 



The name "generalized electroacoustic reciprocity lemma 1, seems appropriate for (35) since 
it represents a merger of the Lorentz lemma [16| of em theory and the acoustical reciprocity 
theorem ((45) below) as applied to adjoint systems. For self-adjoint transducers, i.e , transducers 
for which the transposed material parameters of (30a) equal those of (30b), (35) reduces to the 
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corresponding Primakoff-Foldy expression [7] in which a second set of fields on the same transducer 
replace the adjoint fields. Auld [12 1 uses this latter expression with the complex conjugate fields to 
derive the orthogonality relation [18J between different modes of piezoelectric waveguides and 
resonators. Equation (35) demonstrates that the waveguides and resonators need not be restricted 
to piezoelectric material in order for the modes to satisfy orthogonality. 

A more general expression similar to (35) may be derived from the equations represented 
by the operator L. Each row of the operator L in (24) could be multiplied by an arbitrary constant 
without affecting the essential content of the equations. Specifically, if the first, second, and third 
row of L are multiplied by arbitrary complex constants, Ci, c 2 , and c 3 respectively, (35) is replaced 
by 



I [d(Ex //«)- c 2 (E a X H)]-nda + c 3 | (T a -u-T'a a )'hda=0. 

JAq J A 



(36) 



Although (36) has mathematical validity and interest, the constants d, C2, and c 3 must be 
restricted if (36) is to be used in the derivation of the generalized reciprocity relations between a 
transducer and its adjoint transducer. (The adjoint fields will depend upon C\ 9 c 2 , c 3 .) First, be- 
cause the right side of (36) is zero, the constant ci, can be set equal to unity without loss of gener- 
ality. Second, c 2 must equal C\ = \ to accomplish the derivation of the scattering-matrix reciprocity 
relations ((53) below). And third, if the adjoint operator underlying (36) is required to describe 
the fields of a second transducer (the adjoint transducer), comparison of (29) with (24) (when the 
constants Ci, Ci, and c 3 are included) shows that for C\ = c 2 = 1, the only allowable values of c 3 
are ± 1 unless both Bo and p^are zero. If both Z? and p% vanish throughout the transducers, the 
adjoint operator represents an adjoint transducer for all values of c 3 . 

Equation (35) already represents the case of c 3 = -f- 1, where the relationship between the 
adjoint transducer parameters and those of the original transducer is given by (30). 

Forc 3 = — 1, the operator/, in (24) changes to L , the superscript ( — ) indicating that the third 
row of L has been multiplied by a minus sign. The operator adjoint to L describes an adjoint 
transducer in which the transposed parameters 

a,, ?/, fit, — ft, — v t , — a f , fit, — a,',/2,',?f,fi ,— p° e , and p° tm (37) 

respectively, replace the original parameters listed in (30b). The transducer of (30a) (c 3 = 1 ) converts 
to the transducer of (37) (c 3 = — 1), and vice versa, by reversing the signs of the off-diagonal 
elements in the third row and column of L" . 

When the transposed parameters (37) assume the values of the original parameters (30b), 
the operator L and its adjoint L " become identical to form a self-adjoint operator for which 
the adjoint and original transducer are one and the same. In that case,p" must be zero (p ( j = — p () ) , 
which corresponds to the magnetic-type coupling of Primakoff and Foldy [7|. 

The equations represented by the operators L ana /> " describe tlie same transducer. How- 
ever, the equations represented by the adjoints L" and L " refer, in general, to different adjoint 
transducers. Section 5 shows that a transducer self-adjoint with respect to L satisfies reciprocity, 
while a transducer self-adjoint with respect to L satisfies antireciprocity. 

5. Generalized Reciprocity Relations 

The derivation of the desired reciprocity relations can now be accomplished with the help 
of the reciprocity lemma (36). It will also be demonstrated that the SIM relations (9) and PWSM 
eqs (8) result as a by-product of the reciprocity derivation. Such a demonstration is included 
simply to prove that the two matrix descriptions of a transducer, which satisfies eqs (12-19) , need 
not be postulated as a definition of linearity, but follow from the linearity of (12-19). 
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The reciprocity lemma (36) for adjoint transducers may be written as 

I [(E"XH) — (ExH a )]-hda± ip a u n — pu«) da = 0, (38) 

Ja Jaq 

where the plus sign refers to the operator /. and the minus sign to L . The surface Ao of the trans- 
ducer has been taken in the surrounding fluid, so that the stress tensor relates to pressure like 

— h - T - h = p 

— h - T" • h = p" , 

except on the feed area .So where the contribution to the integral is assumed negligible. 
The integral overdo of the em fields expands to 

Vff* - Voli + I [(E" x H) - (Ex H")] • hda (39) 



f [(E"xH 

J Aq-Sq 



(h points into the transducer), after substitutions of the modal fields normalized to unity by (7). 
The contribution from the integral over Aq — So, the area of the transducer excluding the area 
So, vanishes if either of the following conditions is satisfied: 

(1 ) The surface A {) — So of the transducer (but not necessarily the surface A p — So of the power 
supply) is electromagneticall) shielded so that E. E" or//, H" are zero on Aq — So. (Prescribed 
boundary conditions on t he adjoint transducer are chosen to follow those of the original transducer. ) 

(2) The surface A p — So of the power supply (but not necessarily the surface Aq — So of the 
transducer) is electromagnetically shielded and no em sources exist external to the transducer- 
plus-(power supply or detector) A. (Regions of linear em material are permitted outside A,) Under 
these conditions, the surface A {) — So may be enlarged to a surface of infinite radius where the 
Sommerfeld radiation condition or the exponential decay in electromagnetically lossy media de- 
mands that the integral be zero. 

11 the transducer operates at quasi-static frequencies so that the em fields outside the trans- 
ducer obey 

£ = -Vi//,£'' = -Vi//'' (40a) 

vxfl=;,vxff«=j", (40b) 

modal theory may not apply, but substitution of (40) into the em integral of (38) shows that the 
contribution overdo - So still vanishes. That is, no shielding is required at quasi-static frequencies 
of operation for the expression 

V'ih-Vol'i 
to represent the em term of (38). 

With the em integral overdo— So zero, (38) becomes 



V'£lo-VolZ± \ (p a u n -pal) da = 0. (41 

J <<0 



Because the fields of the transducer and its adjoint exist ii lependently and the velocity-current 
variables of each transducer may assume arbitrary values, (41) transforms into the SIM relations 
by first choosing 

/;;, = 1, «;;,(r) = o, (42a) 

which gives 

V = V'£Jq ± p ( {(r)u n (r) da; (43a) 

J, o 
then choosing 

/;;., = (), //',;,(f)=o(r-r'), (42b) 
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which gives 

p(r) =± V; Y> (f)/ H- p'l (f,f )tt„(f«) da<>. (43b) 

J.lo 

Comparison of (43) with (9) shows that the various spatial impedances equate as 

Zb=Vl x , h' (f) = ± p« (?) , 

(44) 
A(f)=±r& ( f )' Zo(f,fo)=pg (r,r () ). 

The PWSM equations (8) are disclosed in a similar manner once (2, 3, 6) and the acoustical 
reciprocity theorem, 

J (p'ttJ-pX) da + j ( P l iK-p-u{) dR = 0, (45) 

as applied to any two acoustic fields supported by the fluid between A and Si , convert (41) to 

iyo(6oag-6ga )± f r)(K) [a(K)b« (-K) -a" (-K)b(K)] dK=0. (46) 

K 

Since the values of a , a and a{{, a" may be designated arbitrarily and independently of each 
other, choose 

a'i,= l, a'J(£)=0, (47a) 

and 

a ( '; 2 = 0,a 2 '(£) =S(JCH-Z), (47b) 

to extract the PWSM equations from (46), 



6o=i;; i ao+- r){L)b ( {{-L)a{L)dl (48a 

7 

6(X)= + -^-6;'(X)a(, + -^- f v (L)b"(K-L)a(L)dL. 

17(A) ~ 17(A) J/. 



(48b) 



Equations (48) coincide with (8) under the replacements, 

Soo=6;; 1 ,So,(L)==hr;(L)6; , (-L)/T;o, 



(49) 



If the conditions of (42) and (47) are applied to the given transducer rather than its adjoint 
transducers, (41) and (46) produce (43) and (48) with the given transducer and adjoint transducers 
interchanged. That is. the SI M relations and the PWSM equations describe the adjoint transducers 
as well as the original transducer. 

The generalized reciprocity relations between a given transducer and its adjoint transducers 
may also be extracted from (41) and (46) by eliminating the pressures and voltages of (41) with 
the SIM relations, and the outgoing amplitudes of (46) with the PWSM equations. Such a procedure 
yields 
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/J!(Z;;-Z/, )/o + /o [A w '(f)+A(f)]a;j(f)da-/g [A'(r) + A w (f)]Mr)da 

± I «;;(f)[Z;;(fo,f)-Zo(f,f»)|w w (f«,)da»da = (50) 

J lu J I,) 

for (41), and 

Tjoag(Sgo-Soo)ao + a« J [ijoSg, (L) =F7 ? (L)Sio(-L)]a w (L)dZ 

+ «S f [7?oSoi(Z)+7?(L)Sf (-L)]a(Z)dL (51) 

± f [ a"(L)[r ] (K)S ( ; i (-K,L)-r ] (L)Sn(-L,K)]a(K)dKdL = 

for (46). 

Since they must hold lor arbitrary values of the inputs, (50) and (51 ) imply 

Z'; = Z h (52a) 

h"'{r)=±h{r) (52b) 

A"(f)=±A'(f) (52c) 

Zg(f , f)=Z (f, f ) (52d) 

and 

Sq == Soo (53a) 

i)oSg 1 (K)=±i|(K)S,o(-K) ( 53b ) 

i|(J5:)Sf (K)=±ijoSoi(- J K) (53c) 

i){K)S'l x (-£ L)=r 1 (L)Su (-1, K) . (53d) 

Either tin 1 "6" or "c" equation of (52) and (53) may be considered redundant since the entire 
derivation could be repeated with the given and adjoint transducers interchanged. 

It is (53b) that enables the characteristic receiving matrix Soi of an arbitrary linear transducer 
to be expressed mathematically as a function of the equivalent sources within the adjoint trans- 
ducers [19|. As mentioned in the introduction, such expressions have special utility in the extrapola- 
tion techniques of the PWSM formulation [4, 8]. 

When a transducer is self-adjoint with respect to L or L~, (52) and (53) reduce to the original 
reciprocity relations of (10) and (11). Reciprocity defined by both Kerns and Foldy-Primakoff is 
satisfied by an electroacoustic transducer that exhibits electric-type coupling (e.g.. the condenser 
microphone) to form a self-adjoint operator />. Antireciprocity is satisfied by a transducer that 
exhibits magnetic-type coupling (e.g.. the dynamic speaker) to form a self-adjoint operator L . 
If the static fields as well as the piezoelectric and piezomagnetic coupling coefficients vanish, 
comparison of (30a) with (37) reveals that both adjoint transducers become identical to the original 
transducer, and by (52) and (53), the transmitting and receiving parameters h, h' , Sio, Soi also 
vanish. In other words, there would be no electroacoustic coupling— a result that the constitutive 
relations (14-18) confirm. 
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6. Adjoint Reciprocity Theorems 

The generalized reciprocity relations (52) and (53) may find expression in a variety of derivative 
"reciprocity theorems and principles." As a consequence of (52) or (53), the "electroacoustic 
reciprocity theorem" [6, 20 1 for reciprocal transducers extends to nonreciprocal transducers by 
the relating of the microphone response of the given transducer to the speaker response of its 
adjoint transducers (and vice versa). Similarly, the "principle of reciprocity" [21] that relates 
scattered pressure at a point B from a point source at A to scattered pressure at points when the 
point source is moved to B may be extended to nonreciprocal transducers through the concept of 
the adjoint transducer. The following derivation of the aforementioned results are performed 
easily and in great generality by using the Kerns PWSM description. 

The standard speaker response S(r) of a given transducer is defined as the ratio of the pressure 
at a point r in the ambient fluid to the transducer input current /<> [6], 

C(-\ -Rill ,rA^ 

S(r) = —T-. (54) 

Because incoming acoustic waves are assumed zero, i.e., a(K) = 0, the pressure in (2a) may be 
expressed in terms of b(K) = V(K) alone, 

p(r)=7T \ b(K)etv*eX''dK. (55) 

The factor e iyz appearing in the integral of (55) allows the z-coordinate of Fto assume values other 
than zero. (The fluid must contain no sources or inhomogeneities between z=0 and the z in (55).) 
The PWSM equations (8), which reduce to 

&o = Sooa<) (56a) 

b(K)=Sio(K)a« (56b) 

when a(/C) = 0, combine with (6b) to yield the following expression for b(K) in terms of/ : 

b(K)= f f ,,,( f(" ■ (57) 

(1 — OOOJTJO 

After substitution of b(K) from (57) into (55), the speaker response (54) takes the final form of 

1 



S(r) = - - \ ■ ) S,»(K)e'*' ? dK,k = K + ye z . (58) 



The adjoint microphone response is defined as the ratio of the open-circuit voltage Vo for an 
adjoint transducer in the presence of a spherical pressure wave centered at point f in the fluid to 
the pressure in the spherical wave at a reference point f c [6], 



Pr(rc) 


(59) 


j/e*|f c -f| 


(60) 


Pr\r,-)— . _, - 
\Tc-T\ 
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(The amplitude of the spherical wave is denoted by j#.) Since /<> is zero, a () equals bo, and Vo is 
equal to 2ao. The scattering-matrix equation (8a) for an adjoint transducer relates a<> to a(K) , 



<>»^ {l -! S: ,j JK(K)a(K)dK. 
The incoming amplitudes (i(K) are given by (4a) as 

a(K) = F(K) =^- f p r (« r )e-^«' rf/? r . 

277 J.v, 



(61) 



(62) 



Replacement of the pressure in (62) by its value given in (60) enables the explicit determination 
ofa(K), 



Mi 
a(K) = — e iyz e-*R, 

y 



(63) 



which combines with (61), (60), and (59) to produce the final form of the adjoint microphone 
response, 



lid [ s;; 



(-K) 



y 

Since the adjoint reciprocity relations require that 

S'L — So 



eikrdK,d= \fr-r\. 



(64) 



T/oS^(-^)=^^S 1 o(^), 
COpo 



(58) and (64) show that 



M"(r) 
S(r) 



(Opo 



(65) 



or that the magnitude of the ratio of microphone response at r of the adjoint transducers to the 
speaker response at r of the original transducer is a constant independent of the position r and of 
the particular transducer involved. The constant could also be made independent of the distance 
d to the reference point f c if the microphone response were 'normalized" to the amplitude ,s/ of 
the spherical wave. If the transducer is self-adjoint (reciprocal or antireciprocal) (65) becomes 
identical to the standard electroacoustic reciprocity theorem [20 1. 

Next, consider the derivation of the adjoint "principle of reciprocity" which applies to the 
scattering of a spherical wave from a given scatterer (or transducer terminated in a passive load). 
The scattered pressure p S c{fn, F.\ ) at a point r H in the fluid caused by a spherical wave centered 
at a point Fa in the fluid is given, as in (55). by 



p„.(f„,r,) = J- f b(K)e' kr dK, 

277 J, 



(66) 
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where b (K) may be found in terms of a (K) from the scattering equations (8): 

b(K)= l L { S '°fl ) -7ls!!) r +Sn (K:L))a(L)dL. (67) 

The reflection coefficient Ft is related to the terminating load Zt = ~~ 'Vol 7o by 

Vr -f^. (68) 

Since the source is a spherical pressure wave centered at r A , the amplitudes a(L) take the same 
form as in (63), i.e., 



a(L)=—-r-e e 69 

Substitution of (69) and (67) into (66) results in the scattered pressure 

/ S to {K)S ot {-L)%) | , i|(/? _ Z) \ 
Psc (n , ri ) = — ^ -j- e lk rB e d rA dL dK. (70) 

Z7r J k jl y\L) 

If the position of the source and observer is interchanged and the given transducer is replaced 
by an adjoint transducer, the scattered pressure at T\ is found in a similar manner to be 



suK)s%(-L)r T . 



+ S« 1 (K,-L)) 

' ik ■ 



(71) 



p" (r A9 r B ) = — f f - e m ' rA e u ' r ' dL dK. 

2irJ*Jl y(L) 

(The termination for the adjoint transducers is kept atZr. ) 

The adjoint reciprocity relations (53) demonstrate from (70) and (71) that 

Psc(f B , r A ) = p a sc (r A , r B ) , (72) 

which completes the proof of the adjoint principle of reciprocity. (Note that the bracketed quantity 
in the integrands of (70) and (71) represents an equivalent scattering matrix with reciprocal proper- 
ties identical to Sn.) That is, the scattered pressure measured at a point r B from a point source at 
?a is equal to the scattered pressure measured at the point ?a when the point source is moved to 
r B and an adjoint transducer replaces the given transducer. If the transducer is self-adjoint (re- 
ciprocal or antireciprocal) (72) expresses the usual principle of reciprocity [21] for scatterers. 

7. Power Relations 

In this section power relations are developed and used to show that the dissipative character- 
istics (lossiness, losslessness, and "gaininess") of the material of the given transducer are preserved 
pointwise in the adjoint transducers. Such a result lends support to the feasibility of producing the 
adjoint transducers in the laboratory. 
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For harmonic fields, the time-average em and acoustic power flow per unit area is given by 
-Re [E X //*] and — — Re [T • /7 :f: ] respectively. (The asterisk denotes the complex conjugate.) 

The total time-average power/ 3 flowing into the electroacoustic transducer through its surface Aq 

is, therefore, given by 

P=oRe<j> [E x H [ -f •//*] >nda, (73) 



where ft denotes the unit normal into the transducer. Since the surfaced lies in or on the boundary 
of the ambient fluid, 

n • T - //* = — pu n *, 
and the power equation (73) becomes 

P = |Re<i> [(EXH*) -n + puf] da. (74) 

If the transducer is electromagnetically shielded, the em part of the integral (74) evaluates simply 
as VJ* by expanding^ and//* in terms of the waveguide mode, i.e., from (5) 

Et=V e (h H? = Io*ho, 

and making use of the normalization equation (7). (Waveguide theory shows that the basis holds 
(e<), hi) ) of a propagating mode in a uniform, isotropic guide may be chosen real.) 

The SIM relations (9) convert the power expression (74) to an Hermitian form containing con- 
tinuous (u,i(T)) and discrete (Jo) variables: 

P = /„(Z„ + Z*) + f [h{h'*(r) + h{r))u*(r) + I*(h' (f) + h*(r))u„(r)] da 

f f (75) 

-f tt„(fo)[Zo(f, fo) +Z*(f , r)]u*(r) da {) da. 

JAq JAq 

(The taking of the real part has been accomplished by adding each term to its complex conjugate.) 
A transducer may be classified as lossy, lossless, or "gainy" if, for all values of u„ and / () 
(except u„, h) all zero), the power P is positive, zero, or negative, i.e., if the Hermitian form (75) 
is positive definite, zero, or negative definite [22]. 

Mercer's theorem [23], which can be generalized to provide necessary and sufficient conditions 
on the spacial impedances Z&, h, h' , and Z» for the power P to belong to one of the "definite" value 
classes, reduces to very simple relationships among the spacial impedances when the transducer 
is lossless. In the lossless case, i.e., when the power expression (75) is identically zero for all values 
of u n and/o, various combinations of / () equal to zero and u„ (f) equal to zero or a sum of delta func- 
tions can be chosen to prove that the spatial impedances satisfy the following equalities: 

Zb + Z * = (i.e., Z b is imaginary) (76a) 

h'(r)+h'Hr) = (76b) 

Zo(f, r )+Z*(fo, f) = 0. (76c) 
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A reciprocal or antireciprocal transducer also satisfies (11), which combine with (76) to show 
that Zi\ (as well as Zi,) is imaginary, and that both A, h' are imaginary or real depending upon 
whether the transducer is reciprocal or antireciprocal respectively. It can be demonstrated [9] 
that no simple relations similar to (76) exist among the scattering matrices of lossless transducers. 

Substitution of the adjoint relations (52) for the spacial impedances into (75) yields an expres- 
sion (similar to (75)) which can be manipulated slightly to prove that the power input to both adjoint 
transducers (one associated with L, the other with L ) belong to the same value class, i.e., have the 
same definite or indefinite form, as the power input to the original transducer. However, the 
power input to the adjoint transducers does not, necessarily remain the same as the power input 
to the original transducer for given excitation /<» and u,,(r). 

Finally, consider the problem of expressing the power relation (73) in terms of the internal fields 
and material parameters of the electroacoustic transducer. The divergence theorem transforms the 
surface integral (73) into a volume integral over the volume V o (enclosed by So) of the transducer 
to produce, 

P=\ [H* • (Vx£) + ^(V XZ?*) -l-(Vx #*) -E- (V x H) 

(77) 

- (V ■?) • a*- (V •?*) • il-f : V/!*-?* : Vfi] dV, 

after the divergence is expanded, and the taking of the real part is accomplished by adding each 
term to its complex conjugate. The equations (12-19) allows the integrand P v of the power integral 
(77) to be written as an Hermitian form, 

[//*, £*, s*, a*]/ia)(jjL-JjL*) i(o{v-T*) M/J-ft'*) 

ico(f — v*) [io)(i — i*) — (o"+o"^*)]ico(a — a' t ) BnX (a+d - * 

Mj3'-j3*) iv(a'-a*) -w(c-c*) 

-B X (cr+5,*) B„X (a+aflXft, 



Pr- 




or simply 



P r = X*%TX. (79) 



P r may be interpreted as power input per unit volume of the transducer. 

A sufficient condition for the total time-average power input to be positive definite, zero, 
or negative definite is for the Hermitian integrand (78) to be correspondingly positive definite, zero, 
or negative definite. Likewise, a standard theorem of linear algebra [22] states that the matrix ffl 
must belong to the same value class as the associated Hermitian form (78), provided //, E, 
s and u may assume arbitrary values at a point. In particular, if P v is zero, the transducer is 
lossless, cW must be zero and therefore the material parameters satisfy the following simple 
relationships: 

cr = —crf. € = €*, /1= /xf, i> = ff, c = cf, a' = af, and /§ ' =/3 r * (80) 

If, in addition, the transducer is reciprocal or antireciprocal, (30) and (37) show that the 
material parameters also satisfy 

o-=cr/. e = e f . Ji = Ji t , ~v = Tt, c = c t , a' =±a t , and f3' =+/3,* , (81) 

where the upper and lower sign in the last two equations of (81) refer to reciprocal and antireciprocal 
transducers respectively. Equations (80) and (81) together imply that a lossless (in the sense of (80)) 
reciprocal or antireciprocal transducer possesses real o~, 6, jlx, c, and imaginary v and r. The 
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piezoelectric triadics a, a' are real for reciprocal and imaginary for antireciprocal lossless trans- 
ducers. Conversely, the piezomagnetie triadics /3,/3' are imaginary for reciprocal and real for 
antireciprocal lossless transducers. 

Equations (30) and (37) may also be used to demonstrate from (78) and (77) that the total power 
input for the two transducers adjoint to the original transducer belong to the same value class or 
have the same dissipative character as the original transducer — a result which concurs with the 
previous analysis on the spacial impedance power relation (75). 

8. Concluding Remarks 

A brief summary of the main properties of the L and L adjoint transducers is displayed in 
tables 1-4. From an experimental standpoint, it is naturally of interest to know if the adjoint trans- 
ducers are physically realizable, i.e., producible in the laboratory. For example, if a transducer 
adjoint to a given transducer could be produced, the receiving characteristics (Soi) of the given 
transducer could be determined by measuring the transmitting characteristics {S%) of the adjoint 
transducer. Such an indirect experimental procedure would be desirable with a set-up designed 
to measure transducer far-field patterns. 

Table 1. 
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€ 
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Mass density 




Pi 
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Static magnetic induction 




So 


k = - b° 




B''= Bo 


Static charge distribution 




Pi 


P? = Pe 




Pe° = ~ Pe 



Table 2 



(oven transducer 
(described by eqs (12-19)) 


L Adjoint transducer 


/. Adjoint transducer 


External characteristics 


External characteristics 


External characteristics 


Transmitting (IWSM) 

Receiving (PWSM) 
Acoustic scattering (PWSM) 

EM reflection (PWSM) 


Sio(K) 

Soi(K) 
Sn(K,L) 

Son 


SM/C)=So,(-/Ot,o/t, 

S[' n (K)=S ] ,(-K) v h } o 
S« n (K,L) = 

Sii(-L t -K)7i(L)lri{K) 


S« lQ (K)=-Soi(-R)7i<J-n 

S'^(K) - = S w (-K)r } h)o 
S t i(-L,-K)T,(L)lr)(K) 

Sqq = S<><) 
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Table 2 -Continued 



Given transducer 
(described byeqs (12-19)) 


L Adjoint transducer 


L Adjoint transducer 


External characteristics 


External characteristics 


External characteristics 


Speaker transfer {SIM) 

Microphone transfer {SIM) 

Open-circuit acoustic {SIM) 
Blocked electrical {SIM) 


h(r) 

h'{r) 
Za(r. f<>) 

z b 


h"(r) = h'(r) 

h"'(r) =h{f) 

Z\\(r. f«) =Zo(f<>, r) 
Z'fl = Z b 


h«(f)=-h'(r) 

h"'{r) =-h{f) 
Zl(r, F ) =Z (fo, r) 

z° = z„ 



Table 3 



Given transducer 
(described byeqs (12-19)) 


L Adjoint transducer 


L Adjoint transducer 


Power input 


Power input 


Power input 


Passive transducers p > 

p ^0 
p = 


II W V 

ooo 


p" > 
p" ^ 
p" = 


Active transducers p < 

/; ^ 
Indefinite 


/;" < 
p" ss 
Indefinite 


p" < 
p" ^ 
Indefinite 



Table 4 



( riven trai 
(described by 


sducer 
jqs (12-19)) 


L Adjoint transducer 


L Adjoint transducer 


Adjoint recipnx 


ity theorems 


Adjoint reciprocity theorems 


Adjoint Reciprocity theorems 


Speaker Response 


Sir) 


S"(r) = M{r)<i>p»e ik ''l47Tid 


S"(r) =-M{r)a)p l) e ikd l47rid 


Microphone response 


M(r) 


M"{r) =S(r)47Tidlwp»e ih " 


M"{r) =-S(r)47ridl(op tt e ikd 


Scattered pressure 
at f H from point 
source at r,\ 


Psc(r~n* r.\) 


P'scirn, r.\) = Ps<(r\* fa) 


P's C (rn, r.\ ) = PsAr.\,r H ) 



The complex problems involved in synthesizing transducer materials of prescribed electro- 
acoustic properties is not the subject of the present paper. However, two groups of nonreciprocal 
transducers should be mentioned for which the creation of the adjoint transducers require, in prin- 
ciple, a relatively simple procedure. Transducers in both groups have symmetric conductivity, (7, 
permittivity ?, permeability /Z, Hooke's tensor c; and magneto-electric tensors related by F=— r*. 
Transducers of the first group have, in addition, the piezoelectric and piezomagnetic tensors 
related by 



■fi\. 



Transducers of the second group satisfy 



Table 1 shows that the L adjoint to transducers in the first group can be realized by simply 
reversing the direction of the static magnetic induction B<> which biases the given transducer. (For 
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some gyrotropic media, e.g., ferrites. transposition of the permittivity and permeability tensors is 
also accomplished by the reversal of the magnetic induction. In fact, the reversal of Bo to form the 
L adjoint transducer is analagous to the reversal of the dc magnetic field necessary to obtain purely 
em reciprocity relationships for gyrotropic media |24|. ) 

Table 1 also shows that the L adjoint to transducers of the second group can be realized simply 
by reversing the sign of the static charge distribution p". or equivalently, the static electric field 
and electric polarization which biases the transducer. 



The author is grateful to D. M. Kerns for the original suggestion of the work and for the many 
helpful discussions throughout its completion. 
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